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CHAPTER  I 
INTRODUCTION 


Adaptive  arrays  have  been  under  study  in  recent  years  as  a 
means  of  protecting  radar  and  communication  systems  from  interference. 
These  arrays  are  based  on  the  original  work  of  Applebaum  and  Widrow  et 
al.  Applebaumfl]  presented  an  array  control  loop  that  maximizes  a 
generalized  signal -to-noise  ratio  (SNR).  Widrow  and  his  co-workers[2] 
presented  the  least  mean  square  (LMS)  error  algorithm,  based  upon  the 
method  of  steepest  descent.  Both  the  Applebaum  array  and  the  LMS  array 
have  found  extensive  applications  in  radar  and  communication  systems. 

One  of  the  problems  in  applying  adaptive  arrays  to  communication 
systems  is  that  the  array  speed  of  response  varies  with  signal 
strengths.  The  array  speed  of  response  is  determined  by  the  eigen¬ 
values  of  the  so-called  covariance  matrix,  which  is  the  matrix  of  the 
cross  products  between  the  array  element  signals.  These  eigenvalues 
depend  on  signal  powers.  A  strong  signal  produces  a  large  eigenvalue 
and  a  weak  signal  produces  a  small  eigenvalue.  If,  for  example,  the 
array  must  null  interference  40  dB  above  thermal  noise,  the  largest 
eigenvalue  will  be  approximately  104  times  larger  than  the  smallest 
one.  It  is  important  to  keep  the  range  of  variation  of  the  eigenvalues 
as  small  as  possible. 

The  eigenvalues  not  only  depend  on  signal  strengths,  they  also 
depend  on  the  signal  arrival  angles,  the  signal  bandwidths,  and  the 
array  parameters  (element  spacings  and  element  patterns). 

The  purpose  of  this  report  is  to  investigate  and  characterize 
the  actual  behavior  of  the  eigenvalues  in  some  simple  adaptive  arrays. 
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We  examine  here  arrays  with  up  to  four  elements  and  determine  the 
exact  eigenvalue  behavior  as  a  function  of  signal  strengths,  band- 
widths,  angles  of  arrival  and  array  parameters. 

Although  the  problem  of  eigenvalue  spread  in  adaptive  arrays  is 
well  known,  very  little  data  exists  in  the  literature  showing  actual 
eigenvalue  behavior.  Some  information  on  eigenvalues  has  been  given 
by  Gabriel [3],  who  discusses  this  subject  in  connection  with  retro- 
directive  eigenvector  beams.  His  paper  also  gives  some  qualitative 
descriptions  of  the  effects  of  eigenvalues  on  array  performance. 
Mayhan[4]  has  also  presented  some  data  on  eigenvalues  for  multiple 
beam  antennas.  He  has  considered  the  eigenvalues  for  non-zero  band¬ 
width  signals  by  regarding  the  bandwidth  as  a  perturbation  of  the 
original  CW  covariance  matrix[5].  However,  these  papers  do  not  give  a 
complete  overview  of  eigenvalue  behavior  as  a  function  of  the  signal 
and  array  parameters.  Our  purpose  here  is  to  provide  such  data  for 
some  simple  arrays. 

We  begin  in  Chapter  II  by  establishing  notation  and  formulating 
the  problem.  We  then  derive  the  eigenvalues  for  the  case  of  two  in¬ 
coming  signals,  one  desired  and  one  interference.  To  simplify  the 
problem,  we  first  work  out  the  solution  for  zero-bandwidth  (CW)  signal 
In  Chapter  II-C,  we  develop  an  interesting  relation  between  the  array 
output  signal -to-interference-plus-noise  ratio  (SINR)  and  the  eigen¬ 
values.  We  show  that  when  the  interference  is  very  strong,  the  array 
output  SINR  is  equal  to  one  of  the  eigenvalues  less  one.  Chapter  III 
presents  numerical  results  illustrating  the  effects  of  signal  para¬ 
meters  (strengths,  arrival  angles  and  bandwidths)  and  array  parameters 
(number  of  elements,  element  spacings  and  patterns)  on  eigenvalue  be¬ 
havior.  Chapter  IV  contains  the  conclusions. 


CHAPTER  II 

FORMULATION  OF  THE  PROBLEM 


A.  Definition  and  Notation 

Consider  an  N-element  adaptive  array  as  shown  in  Figure  1.  The 
N  elements  are  assumed  to  lie  along  a  straight  line  with  spacing  D 

+■  h  ^ 

between  the  element  and  the  first  element.  The  analytic  signal 
x  (t)  from  the  element  is  multiplied  by  a  complex  weight  w  (t) 

36  36 

generated  from  the  optimizing  network.  The  resultant  products  are 
summed  to  produce  the  array  output  signal  s(t).  For  an  LMS  array. 


the  weight  vector 

W  =  [w15w2..--wn]T  (1) 

satisfies  the  first  order  differential  equation 

+  k<J>  W  =  kS  (2) 

where  $  is  the  covariance  matrix  of  the  array, 

*  =  E{X*XT}  (3) 

S  is  the  reference  correlation  vector, 

S  =  E{ X*r(t) )  (4) 

and  k  the  loop  gain.  In  these  equations,  X  is  the  signal  vector 

X  =  [x1(t),x2(t),---,xN(t)]T  (5) 
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DESIRED 

SIGNAL 


Figure  1.  An  N-element  adaptive  array  with  two  incoming  signals. 
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r*(t)  is  the  (complex)  reference  signal  in  the  array,  T  denotes  trans¬ 
pose,  *  complex  conjugate  and  E{ - >  expectation. 


As  will  be  seen  below,  <!>  is  not  singular  as  long  as  the  element 
signals  contain  independent  thermal  noise.  Therefore,  the  inverse  of 
*  is  always  well  defined.  Let  this  inverse  be  denoted  by  <J>  ,  and 

then  from  Equation  (2)  the  steady-state  weight  vector  is 

“st  -  r's  (6) 


where  the  subscript  "st"  denotes  steady-state.  The  complete  time 
response  of  the  weight  vector  is  then 


M 

W(t)  =  l  C 

i'=l  ? 


-kx  t 
e  *  + 


(7) 


where  the  C ' s  are  constant  vectors  depending  on  the  initial  condi¬ 
tions  of  the  weights  at  t=0.  The  x{'s  are  the  distinct  eigenvalues 
of  <t>  and  M  is  the  number  of  distinct  eigenvalues.  It  is  these  eigen¬ 
values  that  control  the  transient  response  of  the  array  weights  and 
that  concern  us  in  this  report. 


We  shall  determine  the  eigenvalues  of  $  under  the  condition  that 
there  are  two  signals  coming  into  the  array.  One  is  the  desired  sig¬ 
nal  arriving  from  angle  0^  and  the  other  interference  from  0..  Both 
angles  are  measured  with  respect  to  broadside,  as  shown  in  Figure  1. 

We  also  assume  each  element  signal  contains  a  thermal  noise  component. 
Thus,  the  analytic  signal  behind  the  element  is  written 

x?(t)  =  cfv(t)  +  ^(t)  +  ne(t)  (8) 

where  cff(t)  and  iv(t)  are  the  received  desired  signal  and  interfer¬ 
ence  on  the  element,  respectively,  and  n p ( t )  is  the  element  noise. 


The  above  equation  suggests  that  we  can  divide  the  signal  vector 
X  into  the  sum  of  three  component  vectors,  i.e., 

"  '  V  "t  *  Xn  <9> 

with 

xd  =  [^(t),d2(t),---,dN(t)]T  (10) 

the  desired  signal  vector, 

xi  =  ['i1(t),i2(t),---,lN(t)]T  (11) 

the  interference  vector,  and 

Xn  =  [n1(t),n'2(t),---,nN(t)]T  (12) 

the  noise  vector. 


We  shall  assume  here  that  the  desired  signal,  the  interference 
and  the  noises  are  zero  mean  Gaussian  random  processes  uncorrelated 
with  each  other.  Thus,  the  signal  vectors  are  statistically  independ¬ 
ent  of  each  other,  i.e., 

E{X*xj}  =  E{X*X^}  =  E{X*xJ}  =  E{X*xJ}  = 

Then  from  Equations  (3)  and  (9),  we  have 

*  =  E{(Xd+X.+Xn)*(Xd+X.+Xn)Tl  (13) 

=  EfX*xIi  +  Ef  X*xT  f  +  EO'*xh  (14) 

d  d  it  n  n 

=  <t>  .+<&  .+^  .  (15) 

d  i  n 
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E{X*xJ) 
n  d 


E(<x(> 


=  0 


Consider  first  the  desired  signal.  It  is  clear  from  Equation 
(15)  that  the  desired  signal  part  of  the  covariance  matrix,  4>d,  is 


*  Elxdxd' 


!>„  ] 


em 


(16) 


t  h  t  h 

where  .  denotes  the  matrix  element  of  at  the  f.  n  row  and  the  nr 
d  d 

column.  From  Equations  (10)  and  (16),  we  know  that 


v,m 


Elc^ttfjt)} 


Because  of  the  interelement  propagation  delays,  we  can  write 


07) 


fyt)  =  ft(od)'i(t-Td  )  08) 

where  TA  denotes  the  interelement  time  delay  between  element  and 

element  1,  c((t)  is  the  desired  signal  waveform  and  f  (e)  is  the  volt- 

t  h  ^ 

age  response  of  the  ^  element  to  a  unit  amplitude  test  signal  arriv¬ 
ing  from  angle  n.  In  Equation  (18)  we  have  assumed  that  the  element 
patterns  are  independent  of  the  signal  frequencies  over  the  bandwidth 
of  the  desired  signal  cl(t).  Thus  from  Equations  (17)  and  (18),  we  have 


=  WJUoJEtfMt-TH  )2((t-T,  )} 


?m 


(19) 


The  time  delays  in  the  above  equations  are  determined  by  the  element 
spacings  and  the  signal  arrival  angles,  i.e.. 


T 


d 


i 


sinod 


(20) 


with  c  the  velocity  of  propagation. 


In  order  to  evaluate  the  expectation  in  Equation  (19),  we  use 
the  following  definition.  Let  the  desired  signal  be  a  stationary 
random  process  with  autocorrelation  function 


7 


Rj(t)  =  E{2(*(t)cf(t+t ) } 

(21) 

=  Ef  2f*  ( t-i  )d(t) } 

(22) 

we  can  rewrite  Equation  (19)  as 

f  ,  =  f*(o.)f  (  ',j)Rj(Tj  .  -T  .  ) 

d.  r  d'  m'  d'  d'  d!  d 

(23) 

Furthermore,  we  assume  that  the  desired  signal  has  a  flat,  band-limited 
power  spectral  density  S^(w)  centered  at  w0  as  shown  in  Figure  2. 

Within  the  band  the  desired  signal  has  power  density  2rP^/(At: , ) , 
with  the  desired  signal  power.  Then  the  autocorrelation  function 
is  given  by 


Ry,) 


—  03 


du> 


(24) 


Auj  .  r 

sin  Td~ 

d  Aui^j  i 


(25) 


Combining  Equations  (23)  and  (25),  we  have 

v  =  sinc[f  »“d(Td,-Td  >1 e 

P.m  L  p.  m  J 


H<VTd  i 

z  m 


(26) 


with  sine  x  = 


sinx 


We  can  simplify  the  above  result  by  noting  that  u>  T.  is  just 

o 

the  phase  shift  between  element  1  and  element  z  at  the  center  frequency 


u  .  Let  us  define  this  phase  shift  to  be 


|>  .  ~  «>>  T  , 

dp  0  dp 


(27) 


Upon  substitution  of  Equation  (20)  into  Equation  (27),  we  get 
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desired  signal,  S^(oj). 


2"D„ 

>d,  =  sinod 

«.  o 


with  xQ  the  wavelength  at  the  center  frequency  uQ. 


(28) 


In  addition. 


AwdTd 


t-1 -H?)  (••'.,) -Sv. 


(29) 


where 


Aci)  j 

B  =  ..J 

Q  ni 


is  the  fractional  bandwidth  of  the  desired  signal. 


(30) 


Substituting  Equation  (29)  into  Equation  (26),  we  finally  have 

r,  -1  -*d  ^ 

♦d  -  f;(8d)fm("d)Pd  sinc  ?  BdU'd  -»d  >  e  ”  '  (3,) 

sim  L  f  m  J 

Similar  results  may  be  derived  for  the  interference.  We  have 


i„(t)  -  f,(o.)i(t-T.  ) 

1  '  9. 


(32) 


where  i(t)  is  the  interference  waveform  and 
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D, 


\"r  sinoi 


(33) 


_th 


From  the  previous  discussion,  we  know  that  the  on  element  of  the 
interference  part  of  the  covariance  matrix,  ^ ,  is  given  by 


♦i  =  f  *  ( 0  i )  ^ni  (°i )E^  i*(t_Ti  )i(t-T.  )l 

fm  f  m 

We  define  the  autocorrelation  function  of  the  interference  as 


(34) 


R  > ( t )  =  E{ i*(t)i (t+ i ) 


=  E(i  (t-t)i(t) 


=  P.j  sinC| ^ 


1  n  1  JV 

2  1 1  6 


(35) 

(36) 

(37) 


where  we  have  assumed  the  interference  also  has  a  flat,  band-limited 
power  spectral  density  of  bandwidth  Aw.. ,  as  shown  in  Figure  3.  Hence, 
Equation  (34)  becomes 

♦  ,  =  sinc  [] 


j(<t>-  -4>i  ) 
9.  m 


(38) 


where 


Aw  • 

B-  =  -- 
1  “o 


is  the  fractional  bandwidth  of  the  interference  and 

2"De 

♦i  =  'I—  sinni 
f  o 

,th 


(39) 


(40) 


is  the  interelement  phase  shift  between  the  l  element  and  the  first 
element  for  the  interference. 


Finally,  we  assume  the  noises  are  zero  mean  Gaussian  random 

2 

processes  uncorrelated  with  each  other,  each  with  power  <>  .  Thus, 
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Figure  3.  Power  spectral  density  of  the 
interference,  S-vU). 


<41> 

with  6  the  Kronecker  delta  function.  Therefore,  the  noise  part  of 
the  covariance  matrix  is 

4>n  =  o2I  (42) 

with  I  denoting  the  identity  matrix. 


From  Equations  (15),  (31),  (38)  and  (42)  we  conclude  that  the 


p,mu  element  of  *  is 


c  m 


4>  =  f*(o.)f  (o.)P.  sine  i  B  ,(<j) .  -4).  )  e 

«n  ev  d'  m  d  d  |_2  °  “m  J 

*  Tl  -i 

*  f1(of>fra(»i,Pi  s1"c  [z  e  ' 


j(<t>,-  ) 

m 


+  o2<S 


v.m 


(43) 


In  a  later  section  of  this  report,  we  shall  investigate  the 
eigenvalues  of  <i>  for  arbitrary  bandwidths.  First,  however,  we  shall 
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consider  the  special  case  of  zero  bandwidth  (CW)  signals  because  in 
this  case  we  can  obtain  the  eigenvalues  in  simple  analytical  form. 


B .  Eigenvalues  for  Zero  Band wid th  Si gna 1 s 


For  Bd=Bi.=0,  the  sine  function  in  Equation  (43)  is  unity,  so 
Equation  (43)  simplifies  to 

l(-.i  -+<i  >  .  J(*i  -+i  ) 


Pm  d  i  d'  in  d 
cw 


+  a26 


?,m 


(44) 


This  result  is  equivalent  to  the  form 


4*  =  PM*u]  +  P.U*uT  +  o2I 

cw  d  d  d  ill 


(45) 


where 


Ud  = 


'Jl'd  ‘J  l’d0  ‘Ji'd  “l  T 

fl(,’d>e  •  f.(e.)e  2,  f„(o>  Nl 


N 


and 


Ui  = 


Me^e  f2(e.)e  •••,  fN(e 

are  vectors  that  contain  the  element  patterns  and  interelement  phases 


-J*i“l 

>i  )e  " 


(46) 


(47) 


It  is  helpful  to  work  with  dimensionless  quantities,  and  speci¬ 
fically  to  normalize  the  covariance  matrix  with  respect  to  the  noise 
2 

power  (i  .  We  define 


1 

cw 


cw 


2 

<1 


mV 

dad 


,U*UT  +  I 

l  i  l 


(40) 


where 


1? 


o 


=  the  signal -to-noise  ratio  (SNR)  of  the  desired  signal 
and 


O 


=  the  interference-to-noise  ratio  (INR)  of  the  interference. 

We  shall  determine  the  eigenvalues  of  4>'  rather  than  t  .  The  eigen- 

c.  w  2  w 

values  of  <f>  are  equal  to  those  of  #'  times  a  . 
cw  cw 

Because  of  the  form  of  the  covariance  matrix  in  Equation  (48), 
it  is  clear  that  two  of  the  eigenvectors  of  <t>‘  will  lie  in  the  plane 

CW 

formed  by  and  U*.  Hence  we  may  express  two  eigenvectors  (e)  in 
the  form 


e  =  aUj  +  pU* 


(49) 


where  a  and  b  are  constants  to  be  determined  by  the  requirement 


4>'  e  =  xe  (50) 
cw 

and  x  is  the  corresponding  eigenvalue.  We  may  find  a  and  p  by  sub¬ 
stituting  Equation  (49)  into  Equation  (50).  Straightforward  calcu¬ 
lations  show  that 


>M>  ■  <««M> +  *  ;■ f duduT ) 

*  <(yK  *  s  s^*) 


(51) 


Hence  Equation  (49)  will  be  a  legitimate  eigenvector  if  we  choose  a 
and  p  so  that 
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(52) 


r,.u!l£  +  ~  T.uJu*  =  C.uTu*  +  f  r, .uTUj 

d  d  d  addi  Ji  i  i  f5  i  l  d 
Therefore,  from  Equations  (50),  (51)  and  (52),  we  have 


x  =  l  +  r.  .uju*  +  l-  f.  .uju* 

d  d  d  «  d  d  i 


(53) 


or  equivalently 


J  ■ ' +  f-iu(ul +  f  q"Iud* 


(54) 


By  defining 


£=  Y 


(55) 


and  then  transforming  Equation  (42)  into  a  quadratic  equation,  we  get 


vVX  -  ('•f“K-'-duX)v  -  f,uIud  • 0 


(56) 


This  quadratic  equation  is  readily  solved  to  give  two  solutions. 


1 


1 


(a+b) 


(57) 


Y2  =  TT 

2«Ui 


(58) 


where 


T  *  T  * 
a  =  r.Uiu.  -  T.U'U. 

ill  d  d  d 


(59) 


and 


b  =  (a  +4 r  -  f  A 

id 


uTudj2)1/2 


(60) 
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According  to  Equation  (53),  we  then  have  two  eigenvalues. 


ii =  1  *  ';dudud  *■  \  (a,fb> 


»2  * 1  +  r-duJu3 ♦  \  (a-b> 


Since  we  have  an  N  element  array,  4>‘  is  an  N  by  N  matrix  so 

cw 

there  are  N  eigenvalues.  In  addition  to  the  above  two,  there  are  N-2 
additional  eigenvalues.  To  find  these,  we  note  that  if  is  an 
arbitrary  vector  orthogonal  to  both  and  U*,  then 

*cwei  ■  v^<u3ei> ♦  r-X<u!ei>  ■ («> 

since  Uje-^uje-^O  from  the  orthogonality.  This  result  implies  that  e-j 
is  also  an  eigenvector  of  4’  with  unity  eigenvalue.  In  general,  since 

CW  ^ 

Vw  is  of  order  N,  we  can  find  N-2  such  vectors  orthogonal  to  both  Ud 
and  U*.  Hence  the  remaining  N-2  eigenvalues  are  all  unity. 

We  now  have  found  all  N  eigenvalues  of  4>'  .  In  the  following 

CW 

paragraphs,  we  shall  make  some  observations  about  the  results  obtained. 


First,  all  the  eigenvalues  of  $>'  are  real  because  1  is  a 
a  cw  cw 

Hermitian  matrix.  Moreover, 


A,  >  1  (64) 

for  e=l,2,"*,N  since  both  and  f;.U.U.  in  Equation  (48)  are  non- 

GOG  111  ★T 

negative  definite  matrices.  To  see  this,  consider  for  example  ^d^d* 
It  is  easily  seen  that 


d<u3u>d  -  <'du»d 
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T  -*  * 

where  is  a  non-negative  quantity.  Thus  is  an  eigenvector 
of  fdU*U^  with  eigenvalue  r^uJl^jO.  For  any  other  vector  U  perpen¬ 
dicular  to  llj,  we  have 


(Uju>*  =  (uJu*)U*  =  0  U*  (66) 


*  *  T 

so  U  is  also  an  eigenvector  of  with  zero  eigenvalue.  In  N  space, 
there  will  have  N-l  vectors  U*  perpendicular  to  and  to  each  other, 
and  hence  N-l  zero  eigenvalues.  Therefore  r  u ^llT  is  of  rank  one  and 

T  *  dud-*- 

the  only  non-zero  eigenvalue  is  6^1^.  Thus  is  non-negative 

definite.  Likewise,  '^U*u|  is  also  a  non-negative  definite  matrix, 
so  the  sum 


A  .  tdu*uj  *  qu*uT 


(67) 


is  non-negative  definite.  Since  4>^w  in  Equation  (48)  is  the  sum  of 
the  identity  matrix  and  the  non-negative  definite  matrix  A  in  Equation 
(67),  the  eigenvalues  of  <t^w  are  just  the  eigenvalues  of  A  plus  unity. 
Thus 


*(♦'  )  =  a ( I+A)  =  1+a(A)>1  (68) 

cw 

where  *(•)  denotes  'the  eigenvalues  of'. 

Note  that  in  general  the  number  of  eigenvalues  of  Vw  different 
from  unity  is  equal  to  the  number  of  signals  incident  on  the  array. 
When  the  array  receives  no  signals  other  than  the  thermal  noise,  the 
normalized  covariance  matrix  of  Eqution  (48)  is  simply 


In  this  case,  all  the  eigenvalues  are  unity.  If  one  CW  signal,  char¬ 
acterized  by  strength  >.  and  arrival  angle  n,  is  incident  on  the  array, 
then  the  covariance  matrix  becomes 
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(69) 


*'  =  cU*UT  +  I 

cw 

In  this  case,  one  of  the  eigenvalues  of  is  l+r,U^U*  and  the  remain- 

cw 

ing  N-l  are  all  unity.  From  the  earlier  discussion,  it  is  clear  that 
with  two  input  signals  there  are  two  eigenvalues  different  from  unity. 
In  general,  one  may  show  that  with  K  CW  input  signals  (K<N)  there  are 
K  non-unity  eigenvalues  and  the  remaining  N-K  eigenvalues  are  unity.* 

Note  that  when  only  one  signal  is  incident  on  the  array,  the  one 
eigenvalue  different  from  unity  has  a  simple  form.  Using  the  defini¬ 
tion  of  a  typical  signal  vector  as  Equation  (46),  we  see  that 

1  +  CUV  =  1  +  C  l  |f„(0)|2  .  (70) 

f=l  1 

Hence  this  eigenvalue  is  independent  of  the  element  spacings  in  the 
array  but  is  a  function  of  signal  arrival  angle  o.  Moreover,  if  the 
element  patterns  are  chosen  so  that 

N  2 

i  iyo)i2 

r=l 

does  not  vary  with  o  (as,  for  example,  with  isotropic  elements)  then 
this  eigenvalue  is  constant  for  all  e. 


Next,  we  note  that  A 1  and  in  Equations  (61)  and  (62)  depend 
strongly  on  the  signal  powers.  For  example,  suppose  the  interference 
power  is  much  stronger  than  the  desired  signal  power.  I.e.,  we  have 


and 


(71) 


♦This  statement  assumes  that  all  K  incoming  signals  produce  linearly 
independent  signal  vectors. 
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(72) 


then  in  Equations  (59)  and  (60)  we  may  approximate  a  and  b  by 

a  =  r,.u}u*  (73) 

and 


b  =  a  (74) 

In  Equation  (74)  we  have  neglected  the  term  4^.^  IKU^  c  of  Equation 
(60)  because  it  is  small  compared  to  a,  since 


4 r,  .r.  uV,  2 

id  id 

2 


UTU*  2 
i  d 

“K2 


4r  . 
<  d 


Therefore  the  two  eigenvalues  in  Equations  (61)  and  (62)  can  be  written 
approximately  as  follows: 


The  above  two  formulas  for  the  eigenvalues  depend  solely  on  signal 
arrival  angles;  the  element  spacings  have  no  effect.  It  is  clear  that 
the  largest  eigenvalue  x-j  is  essentially  controlled  by  the  large  in¬ 
terference  and  the  smaller  eigenvalue  is  controlled  by  the  weaker 
desired  signal.  All  other  eigenvalues  are  unity  with  no  dependence 
on  the  signal  arrival  angles. 


Better  approximations  for  the  eigenvalues  x^  and  x^  may  be 
obtained  by  using  the  binomial  expansion  to  approximate  b  in  Equation 
(60).  That  is,  we  use 
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to  get 


1 

i 


(Hx)1/2  =  1  +  i  x 


b  = a  (' +  -P  |uiud|2) 


Then  putting  Equation  (77)  into  Equations  (61)  and  (62),  we  get 

12 


'1 


1  +  r.uTu*  + 

iii  a 


(77) 


(78) 


»2  ;  +  :AUl  -  ^ 


(79) 


The  above  approximations  are  more  accurate  than  Equations  (75)  and 
(76). 


for 

para 


We  now  consider  the  eigenvalues  under  two  special  conditions 
The  first  case  is  when  the  two  signal  vectors  are 
lei ,  i.e. ,  when 


U*  =  h  U*  (80) 

with  h  a  complex  constant.  Clearly,  under  Equation  (80)  we  have  from 
Equations  (59)  and  (60)  that 

a„  ■  [|h|2Ct-rd]uJu*  (81) 

and 

b„  =  K2+4lh|  Vd(U5Ud)2]1/2  *  (82) 

Hence  we  have 
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(83) 


T*  ?  T  *  T*  T  * 

=  l+f.  .U'U  .  +  |h|  r, -ll'l) .  =  1+r  .U'U  .  +  r„U  U. 
d  d  d  1  1  i  d  d  d  d  d  i  i  i 


and 


x  =  1 

mi 


(84) 


In  the  above  equations,  the  subscript  indicates  that  a„,  bH,  A, 

and  A„  are  calculated  under  condition  Equation  (80).  From  our 
^11 

earlier  result  in  Equation  (64)  it  is  clear  that  A,  is  the  smallest 

mi 

possible  value  of 


On  the  other  hand,  if  we  have  orthogonal  signal  vectors,  i.e., 
u{l)*  =  0  (85) 


then  from  Equations  (60),  (61)  and  (62)  we  get 


bA  =  aA 


(86) 


'i,  = 1 


(87) 


A 


2i 


1  *  fdu>d 


(88) 


where  the  subscript  '1'  indicates  that  these  quantities  are  obtained 
under  condition  Equation  (85). 


These  special  case  eigenvalues  specify  the  bounds  within  which 
the  eigenvalues  will  vary  with  signal  angles.  To  show  this,  we  make 
the  following  observations .  From  Equation  (60),  it  is  clear  that* 

b-a  (89) 


T  *  T  * 

*Here,  we  assume  that  a  0,  i.e.,  r -U.U- x  .1)  ,U  ,.  The  case  for  a- 0 
can  be  similarly  deduced.  1 
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Therefore,  from  Equation  (62)  we  see  that 


cduJud 


(90) 


Combining  with  Equations  (68),  (84)  and  (88),  we  see  clearly  that 


'  2 

C II  c 


(91) 


This  inequality  gives  the  bounds  within  which  x ^  varies  as  a  function 
of  and  o..  To  get  the  corresponding  bounds  on  x.j  we  proceed  as 
follows.  Rewrite  Equation  (61) 


x,  =  H^uju*  +  \  (b-a) 


(92) 


Thus  it  is  clear  that 


g  -■  u'iuX 


(93) 


which  gives  a  lower  bound  for  ^ .  Making  use  of  the  Schwartz  inequal¬ 
ity 


T  * 

UiUd 


?  T  **  T  * 

i  (uju,)(udTud) 


we  have  from  Equations  (60)  and  (94)  that 

12 


2  2 

b^  =  a+4.r  .r  . 

i  d 


uTud 


d2*4r(f„(uTu*)(uJu*) 


T  *  T  *  9 

(• iuTui+rdudud) 


(94) 


(95) 


i  .e. , 


b 


'•iuK*'duX 


(95) 
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Equation  (61)  can  aqain  be  rewritten  as 


1 


T*  T*  1  T*  T* 

1  +  '-  -  u:u.+r  .  u'd  .  +  '  (b-r.U  U.-r.  .u'u.) 

ill  ddd  2  iii  d  d  d' 


(96) 


Therefore,  from  Equations  (95)  and  (96),  it  is  clear  that 


I  *  t  * 

l  +  . . u : u .  +  '  u'u. 

tit  ddd 


(97) 


Combining  Equations  (93), 

A,  ,  we  have 
1 1 


(97)  and  the  previous  results  of 

1  ii 


and 


Notice  ttiat  the  eigenvalues  do  not  necessarily  attain  these 
bounds  as  the  signal  angles  vary.  An  eigenvalue,  say  will  often 
have  a  maximum  \„  smaller  than  \0  .  These  bounds  are  determined 

L\m  X  U  ^  y  ^ 

simply  from  the  signal  strengths  and  terms  such  as  U  U  which  evolve 
the  element  patterns  jf  (n)|  but  not  the  element  spacinqs.  Typically 
it  is  found  that  when  the  element  spacings  are  larger  than  half-wave¬ 
length  the  extrema  of  the  eigenvalues  coincide  with  the  bounds  for 
isotropic  element  arrays. 


Next,  we  consider  the  case  of  isotropic  elements.  With  f  (a)=l; 
■’  =  1 ,2,  •  •  •  ,N  we  have  then 

T  *  f  * 

U.U.  =  U-U-  -  N  =  the  number  of  elements  in  the  array, 
d  d  i  i 

Thus  the  eigenvalues  in  Equations  (61)  and  (62)  have  a  simpler  form 


laN  d  i-  ?  (a+b) 

(99) 

HN  d  f  J  (a-b) 

(100) 

where 
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a  >  N(tr/d) 


b  "  (a2*4f.jcd|u[u*|2)'/2  (101) 

The  bounds  in  Equations  (91)  and  (98)  also  reduce  to 

1  +  Nri  I  ^  r.  1  +  N^'i+r(j)  .  (102) 

and 

1  <  Ap  <  1  +  Nf_d  •  (103) 

Now  the  bounds  depend  on  the  signal  strengths  and  ^  only.  It  is 
clear  from  Equation  (102)  that  for  very  strong  interference,  such  that 
and  r  .>>1  hold,  the  variation  in  x^  is  very  small  compared  to 
its  magnitude. 


Since  for  isotropic  elements  the  sum  of  A-j  and  A p  equals  a 
constant,  2+N(f,i+f;d) ,  we  know  that  when  Ap  attains  its  minimum  must 
be  at  its  maximum,  and  vice  versa.  We  shall  show  the  explicit  condi¬ 
tions  under  which  the  eigenvalues  attain  their  extrema  with  respect 
to  signal  arrival  angles.  For  example,  let  us  differentiate  Equation 
(100)  with  respect  to  (K  (with  fixed).  The  result  is 


!i2  =  _  1  *b_  =  flfd  JL  ,|T*  2 

39i  "  2  36.  ‘  b  30.  uiud 


Thus  the  extrema  of  Ap  (with  respect  to  o^)  occur  when 


(104) 


(105) 


or  when 
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m 


3 

30 


T  * 
UiUd 


(106) 


We  see,  then,  that  the  orthogonal  condition  in  Equation  (85)  is  also 
the  condition  for  t°  reach  its  extrema  with  respect  to  o ^ .  On  the 
other  hand,  the  parallel  signal  vector  condition  in  Equation  (80) 
satisfies  Equation  (106)  because 


:)0i 


T  * 
U  U  . 
i  d 


=  !h| 


30. 

1 


T  * 
udud 


Therefore  Equation  (80)  also  gives  the  condition  for  to  be  extrema. 
These  results  indicate  that  the  extrema  of  the  eigenvalues  coincide 
with  the  bounds  given  in  Equations  (91)  and  (98)  for  isotropic  ele¬ 
ments. 


We  now  return  to  arbitrary  element  patterns  and  consider  how 

★ 

the  condition  in  Equation  (80)  can  be  met.  From  the  definition  of  Ud 
and  IK  in  Equations  (46)  and  (47),  it  may  be  seen  that  the  condition 
in  Equation  (80)  requires  suitable  (o^.o.)  pairs  that  satisfy  N  simul¬ 
taneous  equations,  i.e., 

j‘"d, 

h  Vnd)e  =  M°i)e  {107) 

for  y,  =  1 ,2,  ■  •  •  ,N.  One  solution  to  this  is  of  course  with  h=l. 

Whether  other  and  <k  exist  for  which  Equation  (107)  is  satisfied 
depends  on  the  element  spacings  and  patterns.  Actually,  Equation 
(107)  is  just  the  condition  for  a  grating  null [6].  If  the  array  does 
not  have  grating  nulls.  Equation  (107)  will  not  be  satisfied  except 
for  (>.=0^,  so  that  >'s  will  not  attain  the  bounds  in  Equations  (91) 
and  (93). 


Tor  isotropic  elements.  Equation  (107)  reduces  to 


e 


(108) 
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with  h=l .  This  equation  is  always  satisfied  for  o.  =  ed  and  the  symmetry 
angle  o^n-r^  since  the  array  in  Figure  1  has  all  elements  along  a 
straight  line. 


Finally,  we  define  the  eigenvalue  spread  S  to  be  the  ratio 
between  ^  and  Since  both  and  x^  are  functions  of  and  o 
we  have 


x,(o  ,,°  ) 
s(o.,n.)  =  —  -r  - — \ 
d  r  ^2  °d’ei 


(109) 


This  spread  is  then  bounded  by  S  .  and  S  ,  i.e., 

J  mm  max 


mm 


max 


(110) 


where,  from  Equations  (84),  (85)  and  (89) 


'1, 


i+quX 


(with 


UiUd 


=  0) 


(HD 


'll 


max  X/ 


=  1+Cduju*  +  r,.u{u*  (with  U%hU*)  .  (112) 


If,  for  example  we  have  a  three  element  array  with  isotropic  element 
patterns  and  r,d=l ,  £..  =  1000,  we  then  have 


S  •  =  —i-  =  750 

mm  4 


S  =  3004 
max 


The  eigenvalue  spread  is  just  the  time  constant  spread  in  the 
array  transient  response,  as  can  be  seen  from  Equation  (7).  If,  for 
example,  an  array  with  three  isotropic  elements  can  accommodate  a 
time  constant  spread  of  3000  then  the  array  would  be  useful  in  a  signal 
environment  with  interference  power  up  to  ^.  =  1000. 
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In  this  section,  we  have  discussed  the  eigenvalues  of  for 

c  w 

the  case  of  CW  signals.  In  the  next  section,  we  point  out  an  interest¬ 
ing  relation  between  the  second  eigenvalue,  ^  and  the  array  output 
signal-to-interference-plus-noise  (SINR)  for  CW  signals.  Then  in 
Chapter  III,  Section  E,  we  consider  the  case  of  non-zero  bandwidth 
signals. 


C .  A  Relation  Between  SINR  and  the 
Eigenvalues  for  CW  Signals 

Here  we  depart  from  the  main  subject  of  this  report,  the  behavior 
of  the  eigenvalues,  to  discuss  an  interesting  relation  between  the 
eigenvalues  and  the  array  output  SINR. 


With  the  steady  state  array  weight  vector  in  Equation  (6)  the 
desired  signal  component  of  the  array  output  is 

sd(t)  =  «JtXd  .  (113) 

The  output  desired  signal  power  is  then  ("o"  denotes  "output") 

Pod  ■  f  E,l“ItXdl2'  •  (,,4) 

Similarly,  the  output  interference  and  noise  powers  are 


Poi  -  l 


(115) 


and 


4 1 2 


(116) 


We  define  the  array  output  signal -to-interference-plus-noise  ratio 
(SINR)  to  be 

p 

SINR  =  p  .  (117) 

oi  on 
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It  has  been  shown  by  Ishide  and  Compton[6]  that  when  there  is  a  desired 
signal  and  one  interference  signal  and  both  are  CW,  the  output  SINR 
from  the  array  may  be  written 


SINR  =  r 


Jd 


|UX 


i 


(118) 


The  above  equation  is  derived  from  the  matrix  inversion  lemma[7]. 


Note  that  the  vector  products  lllll. ,  uT U .  and  uTu.  which  contro1 

d  d  i  1  id 

the  behavior  of  the  eigenvalues  also  appear  in  the  SINR  formula.  When 

T  * 

the  signal  vectors  are  orthogonal,  i.e.,  U..Uj=0,  the  SINR  in  Equation 
(118)  attains  its  maximum  value.  This  maximum  is 


SINR.  =  r,  .ulUj 
1  d  d  d 


(119) 


T  * 

From  Equation  (85)  we  know  that  1)^=0  is  also  the  condition  under 
which  *2  attains  its  maximum.  On  the  other  hand,  when  the  two  signal 
vectors  are  parallel,  as  in  Equation  (80),  we  find 


SINR,,  =  f. 


d 


T  * 
% 


1H(U>* 


*n2 


S-IM2(u>3 


ri 


(120) 


Then  for  the  case  of  strong  interference,  as  given  in  Equation  (71), 
we  have  SINRM=0.  Therefore,  the  condition  in  Equation  (80)  gives  not 
only  minimum  but  also  a  very  small  SINR.  In  conclusion,  both 
and  the  SINR  reach  their  extrema  under  the  same  conditions.  Hence 
there  appears  to  be  a  close  relationship  between  and  the  SINR  for 
the  strong  interference  case. 
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To  find  this  relationship,  we  now  express  the  SINR  in  terms  of 

the  eigenvalues.  To  do  so,  we  first  expand  the  Hermitian  matrix  <fr' 

c  w 

in  a  spectral  decomposition^] .  We  have 


[*'  ] 
L  cw 


-1 


N 

=  l 

f.=l 


1  *  T 

—  e  e 

A  V  9. 


(121) 


where  the  A, 's  are  the  eigenvalues  and  the  e*‘s  are  the  corresponding 
eigenvectors  of  .  Recall  that  ^  is  given  by  Equation  (48). 

Also  its  first  two  eigenvectors  have  been  found  to  be  (see  Section  B, 
Equation  (49)) 


el  =  VV^i*  =  “1  Ud+BlUi  {]72) 

e2  =  «2(Ud+y2U^  =  “2  Vp2Ui  (123) 

★ 

and  the  remaining  N-2  eigenvectors  are  all  orthogonal  to  both  U,  and 
*  a 
.  Also,  the  first  two  eigenvalues  are  given  by  Equations  (61)  and 

(62)  and  the  remaining  N-2  eigenvalues  are  unity.  As  we  shall  see, 

these  unity  eigenvalues  will  not  appear  in  the  SINR  expression  because 

*  * 

their  associated  eigenvectors  are  orthogonal  to  both  Ud  and  IK . 


First  we  calculate  the  steady-state  weight  vector.  We  assume 
that  the  reference  signal  r(t)  is  a  replica  of  the  desired  signal. 
Then  the  reference  correlation  vector  is  just  the  desired  signal 
vector,  i.e.. 


S  =  yU, 


(124) 


with  y  a  proportional  constant.  From  Equation  (6),  we  have 
W 


st 


)  _]  ★ 
[•!'■  1  yU. 

L  cw  '  d 


(125) 
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Substituting  Equation  (121)  for  [V  ]  ,  we  find  that  the  products 

IT  t  *  cw 

between  and  are  zero  and  do  not  contribute  to  the 

result.  Therefore, 


_  -2 

Wst  '  0  Y 


1  *  T  *  1  *  T  * 

aV  e1  ^e1  Ud  +  T'  e2^e2Ud ^ 


(126) 


A  1  I  I  U 

Since  the  desired  siqnal  is  CW,  we  have 
j(%t+<!> .) 

*d  =  ^ » e  ud 

o 

where  (/Tj  0)  is  the  input  desired  signal  power  P^,  u>0  is  the  signal 
frequency  and  ij/.  is  the  desired  signal  phase.  We  also  assume  that  ^ 
is  uniformly  distributed  between  0  and  2tt.  Now  we  can  calculate  the 
steady-state  desired  signal  component  of  the  array  output  from 
Equation  (113).  We  find 


-  Y/rd  J<“ot+*d> 

sd(t)  »  e 
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From  Equation  (114),  we  have  the  output  desired  signal  power 
P 


od 


2o 


2" 


1 

,,T  * 
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T  * 

i 
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del 
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Ude2 
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Similarly  we  have 


j(wnt+>l'i) 

X.  =  /r"  ,i  e  0  U, 

i  i  i 

o 

where  (/fTo)  is  the  input  interference  power  P.  and  is  the  inter¬ 
ference  phase,  also  assumed  to  be  uniformly  distributed  in  [0,2-n]. 
Therefore,  we  have 


Si(t)  = 


1  |„T  * 

2  1 

T  * 

2 

-  II.  e, 

»,  h  i 

+  - 
X2 

1)  .e0 

i  ? 

(129) 


and 
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The  output  noise  power  is  given  by  Equation  (116),  i.e.. 


,  =  IxL 

on  2fl2 
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From  Equations  (117),  (128),  (130)  and  (131)  we  have 
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(132) 


which  gives  the  SINR  in  terms  of  the  eigenvalues  X-j  ♦  x 2  and  the  pro¬ 
jections  of  each  signal  vector  on  the  eigenvectors. 

It  is  easily  seen  from  Equations  (122)  and  (61)  that 


T  * 
Udel 
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and  similarly 
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Ude2 
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-  Fj1  |  rt^  |  (X2-l  ) 
-1 


=  r 


■]  t,li  I  (  Ai“l ) 
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-  E-  |<»2 j  (x2-l ) 
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(135) 

(136) 


Putting  Equations  (133)  through  (136)  into  Equation  (32),  we  get 
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cd 

cd2,oti l2(Ar1  )2  t  ^2la2l2(x2'1)2 

X1  X2 
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The  above  formula  expresses  the  SINR  in  terms  of  the  eigenvalues  A -j 
and  *2  and  the  coefficients  | n j ' s  and  |  8 1  1  s .  We  need  these  |a.|'s 
and  | 8 • | ' s  to  further  simplify  the  expression. 

In  Equations  (122)  and  (123),  we  need  «-|  and  ^  as  well  as  y1 
and  y2  (given  in  Equation  (57)  and  (58))  to  specify  e*  and  e£  com¬ 
pletely.  The  additional  equation  needed  may  be  obtained  by  normal¬ 
izing  e*  and  e^,  i.e.,  by  enforcing 

T  *  T  * 
elel  =  e2e2  =  1  • 


From  Equations  (107)  and  (123),  we  get 


T-*-’  '2  Al_1  +  l^l2  -7-=  1 

r’d  1  i 


elel  "  |al  1  r. 
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Similarly  we  have 


i^i2  |p2i2 
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Then  from  Equations  (53)  through  (58),  we  have 


(138) 


(139) 
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Substituting  Equation  (140)  into  Equation  (138),  we  can  solve  for 
2  2 

| «i |  and  |p. |  .  The  results  are 


(142) 


Ifq  I2  =  I*!  I2!'*!  I2 


(143) 


with  given  in  Equation  (57).  Likewise,  by  substituting  Equation 
(141)  into  Equation  (139)  we  have 


1  (}  ly2|2  V1 
2-1  w +  “ r ; 


(144) 


and 


I  Bo  I  =  |y?!  I  •’to  1 2 


(145) 


With  the  above  results,  we  may  express  the  SINR  in  terms  of  the 
eigenvalues  and  the  signal  vectors.  The  final  result  is 


SINR 


(146) 


where 
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This  general  expression  is  rather  unwieldy.  But,  for  the  case 
of  very  strong  interference,  the  equation  reduces  to  a  simple  form. 
If  Equations  (71)  and  (72)  hold,  i.e.. 


T  * 

r.  ,u  u. 

i  i  i 


ryA 


and 
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■d 


then  from  our  earlier  results  in  Equations  (75)  and  (76),  we  know  that 
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where  the  first  term  in  Equation  (147)  involving  ^  is  negligible  and 

the  first  term  in  Equation  (148)  is  also  neglected  because  it  involves 

-1  -1 
f.  ,  which  is  much  smaller  than  the  second  term,  proportional  to  r,^  . 

Thus, 
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The  above  expression  can  be  further  simplified  by  recognizing  that 


t . 
'\j  'i 


=  1 


U,-'-  />*/ 

It* 

d|l,dud 

under  Equations  (71)  and  (72).  Therefore  we  have 


SINR  =  \2-l 


(152) 
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This  remarkably  simple  result  is  valid  for  arbitrary  element 
spacings,  element  patterns  and  signal  arrival  angles,  but  only  for  CW 
signals.  The  formula  is  interesting  because  it  tells  us  that  when  the 
interference  is  much  stronger  than  the  desired  signal,  x^  controls  not 
only  the  transient  response  of  the  array  but  also  the  steady-state 
SINR  performance. 

Much  recent  work  has  been  directed  at  the  problem  of  choosing 
element  patterns  in  an  adaptive  array.  The  goal  of  this  work  is  to 
find  element  patterns  for  which  the  SINR  does  not  vary  widely  as  the 
signal  arrival  angles  change.  We  note,  however,  that  because  of 
Equation  (152),  choosing  element  patterns  to  minimize  the  SINR  varia¬ 
tion  also  minimizes  the  variation  in  both  x^  and  the  eigenvalue  spread. 

We  now  return  to  the  main  subject  of  this  report  and  discuss  the 
behavior  of  the  eigenvalues  in  typical  situations. 
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CHAPTER  III 

RESULTS  AND  DISCUSSION 


We  have  determined  the  N  eigenvalues  of  the  covariance  matrix 
for  an  N-element  adaptive  array  with  a  CW  desired  signal  and  a  CW 
interference  signal  incident  on  it.  We  have  also  related  the  eigen¬ 
values  to  the  array  output  SINR,  as  given  in  Equations  (146)  and  (152). 

In  this  Chapter,  we  shall  discuss  the  behavior  of  the  eigen¬ 
values  and  relate  their  values  to  the  signal  strengths  (>'  .  and  . ), 
the  number  of  elements,  the  element  spacings,  the  element  patterns 
and  the  signa1  bandwidths. 

To  illustrate  the  behavior  of  the  eigenvalues,  we  first  present 
Figure  4,  which  shows  a  typical  set  of  eigenvalues;  it  shows  all  three 
eigenvalues  ( \-j ,  and  \^)  versus  the  interference  angle  0.  for  a 
three-element  array.  The  desired  siqnal  angle  is  arbitrarily  fixed  at 
0^=45°  and  the  interference  angle  is  varied  between  0°  and  360°.  Both 
signals  are  CW.  The  elements  are  assumed  isotropic  and  a  half-wave¬ 
length  apart  (D^O,  D2=0.5\o,  D3  =  1 . 0\q)  .  The  signal -to-noi se  ratios 

are  f,,  =  l  and  r,  .  =  10. 
d  i 

From  Figure  4,  we  observe  the  following: 

1)  \,,  the  largest  eigenvalue,  is  always  larger  than  31.  This 
minimum  of  ' ^  is  determined  by  l+.f^UllJ^  from  Equation  (93). 

2)  the  middle  eigenvalue,  varies  between  1  and  4,  i.e., 

1  X2  1  +  d^d^d  ^ron’  ^9,ja^10n  (103). 
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EIGENVALUES 

5.0  10.0  15.0  20.0  25.0  30.0  35.0 


Figure  4.  Eigenvalues  vs  e. 

SNR=1 .0  1 

INR=1 0 

BD=0.0  (DESIRED  SIGNAL  BANDWIDTH) 
BI=0.0  (INTERFERENCE  BANDWIDTH) 
THETA  D  =  45  DEGREES 
3-ELEMENT  ARRAY  WITH  HALF¬ 
WAVELENGTH  SPACINGS. 


37 


3)  \3-l. 


4)  x.|  has  a  range  of  variation  equal  to  3.  i.e  , 

X,  x1  =  34-31  =  3. 
max  min 

This  range  is  determined  by  the  weaker  signal  (the  desired 

signal)  and  has  value  r.U^U^  (from  Equation  (98)).  Also 

we  find  that  x,  =  x,  =  34,  x.  =  x,  =31. 

max  '»  min  1  x 

T  * 

5)  x^  also  has  a  range  of  3,  i.e.,  (From  Equation 

(103)). 


6)  The  eigenvalue  spread  S  has  extrema 

Smax  =  1  +  '-il,X  +  'dUdUd  =  34  (Equation  (112)) 


From  1)  and  2),  we  see  that  the  levels  of  the  eigenvalues  are 
determined  by  the  signal  strengths,  namely,  x-j  by  r  and  x^  by 
x^  is  always  constant  (unity).  Points  4)  and  5)  illustrate  how  the 
ranges  of  the  eigenvalues  are  controlled  by  the  weaker  signal,  in 
this  case  the  desired  signal.  The  graph  also  shows  that  the  eigen¬ 
value  variation  is  symmetrical  around  ik=90°  and  279°,  i.e., 

\(90°-x)  =  x(90°+a)  :  \(279n-t)  =  x(270°+.s)  . 

This  symmetry  comes  from  the  fact  that 


U*(oi)  =  U*(180°-'t. ) . 
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With  this  background,  we  now  discuss  the  effect  of  each  of  the 
system  parameters  separately.  We  start  with  the  signal  strengths. 

A.  The  Effect  of  Signal  Strength 

We  first  show  how  eigenvalues  are  affected  by  the  desired  signal 
strength,  r^.  Figure  5  shows  the  three  eigenvalues  with  f^=3  and  with 
all  other  parameters  the  same  as  in  Figure  4.  (Figure  4  was  for  td=l). 

Comparing  the  two  plots,  we  find 

1)  The  minimum  value  of  A^  is  still  31,  i.e.,  the  minimum  value 
of  A.|  does  not  depend  on  r^,  as  can  be  seen  from  Equation 
(93). 

2)  a^  has  minimum  1  and  varies  between  1  and  10.  The  bounds 
on  Ag  are  from  Equation  (103)  again. 

3)  a3  remains  constant  (unity). 

4)  A,  has  a  larger  range  of  variation,  9  in  this  case,  which 

T  * 

is  CjUjUj  (or  Nsd).  Also  we  see  that 

A,  =  A,  =  40,  A,  =  A,  =  31 . 
max  11  min  1 

5)  Ap  also  has  a  larger  range  of  9,  we  also  have 

A,  Ap  =10,  *2  .  X2„  ' 

max  1  mm  *' 

6)  The  eigenvalue  spread  has  extrema 

S  =40  (from  Equation  (112)) 

max 

S  •  =  =  3.1  (from  Equation  (111)) 

mm  10 
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THETfl  I  (DEG.) 

Figure  5.  Eigenvalues  vs.  o. 

SNR=3.0  1 

INR=10. 

BD=0.0 
BI=0. 0 

THETA  D  =  45  DEGREES 
3-ELEMENT  ARRAY  WITH  HALF¬ 
WAVELENGTH  SPACINGS. 
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7)  The  sum  of  the  three  eigenvalues  is  a  constang  equal  to, 
in  this  case,  42. 

From  the  above  observations,  we  see  that  the  desired  signal 

strength  affects  the  range  of  variation  of  A-j  and  The  range 

is  proportional  to  f.^.  In  addition,  the  increase  in  substantially 

reduces  S  (from  7.75  to  3.1).  At  the  same  time,  increasing  r,, 

makes  S  to  increase  (from  34  to  40).  Therefore,  tripling  f,. 
max  '  3  d 

reduces  S  .  more  than  by  a  half  and  only  causes  S  to  increase 
min  J  max 

less  than  20%. 

The  previous  example  was  for  the  case  of  weak  interference 
(£.=10).  We  now  consider  how  Fd  affects  the  eigenvalues  when  f.  is 
much  stronger.  We  let  f^=1000  and  calculate  x^  for  three  different 
C  -s  (Fd=l,  10  and  100).  The  result  is  shown  in  Figure  6.  We  see 
clearly  that  x^  is  little  affected  by  the  increase  in  F>d  for  strong 
interference.  Thus,  the  larger  eigenvalue  x^  can  be  regarded  as 
independent  of  r_d  as  long  as  f./f,^  >.  10  and  f.  >.  1000. 

In  Figure  7,  we  show  X 2  calculated  for  two  f.^  values  with  all 
other  parameters  the  same  as  in  Figure  6.  Clearly,  X^  depends  very 
much  on  Fd-  Acutally,  for  FJd=10  is  just  a  magnified  version  of  the 
case  Cc|=  1  • 

Next,  we  discuss  the  effect  of  interference  strenqth.  First, 
we  show  that  x^ ,  the  largest  eigenvalue,  depends  mainly  on  the  inter¬ 
ference.  In  Figure  8,  we  have  x^  plotted  against  o.  for  three  inter¬ 
ference  strengths,  f^=100,  1000  and  10000  with  all  other  parameters 
the  same  as  in  Figure  4.  It  clearly  shows  how  for  strong  interference 
the  largest  eigenvalue  x.j  is  a  constant,  as  a  function  of  o  . .  To 
show  the  exact  relation  between  x^  and  f.,  we  present  Figure  9,  which 
is  an  enlarged  version  of  x^(o^)  for  f_=100.  The  percentage  change 
of  x i  from  its  minimum  (301)  to  maximum  (304)  is  so  small  that  x^  can 


I GENVRLUE , 


£i  =10000 


“T - 1 - 1 - I - 1 - j - 1 - I - 1 - 1 - I - 1 - 1111 - 1 — -J - 1 - 1 - I - 1 - 1 — 

90.  180.  270.  360 

THETfl  I  (  DEG.) 

gure  8.  x-j  for  strong  interference. 

SNR=1 . ;  INR=1000,  1000,  10000 
BD=0. 000;  BI=0. 000; 

THETA  D=45. DEGREES: 

3-ELEMENT  ARRAY  WITH  HALf- 

WAVCLENGTH  SPACING. 


indeed  be  considered  as  a  constant  proportional  to  as  Figure  8 
indicates. 

We  now  show  how  the  large  eigenvalue  A-j  behaves  with  weak  inter¬ 
ference.  In  Figure  10,  we  show  A-j  for  f.^=0.01,  0.1  and  1  with  the 
other  parameters  the  same  as  in  Figure  4.  As  may  be  seen,  .\-|  is  nearly 

constant  for  f, .=0.01.  For  r,.<<r.,  we  have 
i  id 

\.j  =  l+3r..  (from  Equation  (93)). 

Then  as  f^.  increases,  A-j  also  increases.  When  r.^>>fc|>  we  have 

A-|  -  1+3/;^  (from  Equation  (93)). 

Thus,  depends  on  r.^  when  F, .  <  ^ r  ^  because  A-j  depends  on  the  stronger 
of  the  two  incoming  signals.  Once  becomes  larger  than  a^  is 
controlled  by  r. . 

Now  consider  how  A^  behaves  as  f._.  is  varied.  Figure  11  shows 
calculated  for  various  interference  powers  from  f. ^  =0 . 01  to  100  for 
the  same  array.  We  first  note  that  A^  approaches  an  upper  bound  as 
\j  increases.  This  upper  bound  is  the  curve  of  output  SINR  less  one 
(see  Equation  (152)).  In  the  figure,  the  top  curve  shows  >2  For 
.  100  and  also  gives  the  output  SINR  from  the  array  (read  from  the 
right-hand  side  scale).  Note  that  the  curve  for  f  .=10  deviates  little 
from  the  •^•100  curve,  so  we  conclude  that  the  approximation 

=  SINR+1 

causes  little  error  as  long  as 

When  .0.01  or  less,  the  second  eigenvalue  is  essentially  unity, 
as  shown.  In  this  case,  the  array  behaves  as  if  there  is  only  one 
signal  present,  the  desired  signal.  Hence  there  is  only  one  eigenvalue 
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Figure  10.  for  weak  interference. 

SNR=1 . ;  INR=0.01 ,  0.1,  1 

BD= 0.000;  BI=0.000; 

THETA  D=  45. DEGREES. 

3-ELEMENT  ARRAY  WITH  HALF- 

WAVELENGTH  SPACING. 
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SECOND  EIGENVALUE  , 

.0.00  1.00  2.00  3.00  1.00 
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Figure  11.  for  different  r . . 
SNR=1 . 0 

INR  IS  CHANGED  FROM 
0.01  TO  10000 
BD=0. 0 
B I =0 . 0 

THETA  D=45  DEGREES 
3-ELEMENT  ARRAY  WITH 
UNIFORM  HALF-WAVE 
LENGTH  SPACING. 
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T  * 

different  from  unity,  \y  The  output  SINR  for  this  case  is  f^d^’ 

from  Equation  (118).  a9  does  not  approximate  1+SINR  because  the  con- 
T  *  i  *  + 

dition  r ,.U .11  - >>f,  .Li  ,U  .  is  not  met.  Thus  the  curves  in  Figure  10  do 
l  i  l  d  d  d  3 

not  indicate  the  SINR  behavior  for  cases  where  r..*  1. 

l- 


From  Equation  (109),  we  see  that  the  eigenvalue  spread  increases 
with  increased  f..  because  is  bounded  by  1+SINR  and  A,  is  proportion¬ 
al  to  r,. .  The  extrema  of  the  eigenvalue  spread  S  are 

1+Nr^  ^  Nr.. 

Smin  ‘  T+NF7  ”  T+nF7 

d  d 


S  ,  =  1+Nr.  .+Nf, , 

max  i  d 


Nfi 


These  extrema  also  increase  with  increased  f. .  The  problem  of  eigen- 

^  3 

value  spread  arises  when  the  interference  is  so  strong  (f ./r^>l 0) 
that  the  array  cannot  accomodate  the  widely  separated  eigenvalues  (or 
the  widely  separated  time  constants).  When  r..  is  small,  the  eigen¬ 
value  spread  is  small  and  causes  no  problem.  Thus  we  are  interested 
in  the  eigenvalue  behavior  under  strong  interference.  For  large  r . , 
A-j  is  essentially  constant  as  signal  angles  and  vary,  so  we  shall 
not  consider  A-j  any  further.  Moreover,  is  unity  for  CW  signals, 
so  the  only  eigenvalue  that  varies  with  signal  angle  is  A„.  In  the 
following  two  sections,  we  shall  concentrate  on  A,,  and  discuss  the 
effects  of  the  number  of  elements,  signal  bandwidth,  etc.,  on  the 
behavior  of  A 


B .  The  Effect  of  the  Number  of  Elements 

In  Figure  12,  we  show  three  plots  of  A^  (as  well  as  SINR)  versus 
o.  with  N=2,  3  and  4.  In  general,  we  obtain  better  output  SINR  with 
more  elements  because  the  larger  N,  the  greater  the  array  gain. 
Correspondingly,  increases  with  increased  N.  Moreover,  as  the 

number  of  elements  is  increased,  the  eigenvalue  spread  also  increases. 
For  example,  the  extrema  of  the  eigenvalue  spread  are: 
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Figure  12.  Relation  between  and  N. 
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UNIFORM  HALT-WAVELENGTH  SPACINGS. 
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=  2003 


S3  =  3004 

max 

S.  =  4005 

Sax 

where  subscripts  2,  3  and  4  indicate  the  number  of  elements.  We 
see  from  Equation  (111)  that  for  isotropic  elements 

r  . 

1  im  s  =  7L 
N-«°  min  d 

T  *  T  * 

because  ILU^U^U^N.  Thus,  the  minimum  eigenvalue  spread  increases 
to  a  limit,  with  increasing  number  of  elements 

Note  that  with  four  elements,  there  are  four  local  minima 
indicated  by  points  C,  D,  E  and  F  in  Figure  12.  For  three  elements, 
there  are  only  two  such  relative  minima  (indicated  by  points  A  and  B). 
Thus,  as  we  increase  the  number  of  elements,  the  number  of  these  mini¬ 
ma  will  also  increase.  The  value  of  at  these  minima  may  decrease 
if  we  change  the  array  parameters,  such  as  the  element  spacings. 

These  minima  may  sometimes  drop  to  unity.  When  this  happens,  there 
is  a  corresponding  null  in  the  SINR.  This  situation  is  due  to  the 

presence  of  a  grating  null  in  the  antenna  pattern[6].  In  addition, 

T  * 

the  maxima  on  these  curves  correspond  to  the  cases  where  IKU^O,  which 
as  we  know  from  Equation  (118),  results  in  maximum  output  SINR. 


S„  =  2?  670 

Sin  J 


S,  =  =  750 

Sin 


S  =  -°-01  =  800 
Sin 


C .  The  Effect  of  Element  Spacing 

Since  the  covariance  matrix  depends  on  the  interelement  phase 
shifts,  which  in  turn  depend  on  the  interelement  spacings,  we  expect 
the  interelement  spacings  to  have  an  important  role  on  the  eigenvalue 
behavior. 
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We  begin  by  considering  a^  for  a  two  element  array  with  several 
values  of  separation  between  elements.  The  elements  are  assumed 
to  be  isotropic.  A  series  of  plots  of  will  be  shown.  The  first  of 

these.  Figure  13.  shows  calculated  for  four  separations,  0,,=0  05\q, 

0 . 1 x Q ,  0.2aq  and  0 . 3 x Q  as  indicated.  In  later  graphs,  is  further 
increased  so  we  are  able  to  see  the  gradual  variation  of  with  D^. 

As  before,  f.^1,  f;.  =  1000,  o.=45°  and  both  signals  are  CW. 


From  Figure  13,  we  see  that  for  very  small  separation  (D^O.05 
A  ) ,  stays  nearly  constant  at  unity.  This  result  indicates  that 
the  output  SINR  of  the  array  is  poor  for  every  n..  Because  the  small 
separation  cannot  provide  enough  phase  shift  between  elements,  the 
array  is  not  able  to  produce  a  satisfactory  SINR. 


o 

o 


o 

O 


SNR  1  ,  INR-1000,  B0-0.000,  BI-0.000, 
THETA  0=45  DEGREES 
? -ELEMENT  ARRAY. 


cr 

z 

H 
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Putting  N=2  in  the  SINR  formula  in  Equation  (118),  we  have  the 
SI NR  for  the  two  element  array  as 


't'i  -4 'A 
2  2  a? 
SINR2  =  2>d  sin  —  ~  ~ 


where  we  have  used 


Ud  = 


- 

- 

1 

U.  = 

1 

■j'I’d 

’  1 

-Jt-j 

2 

2 

Le  J 

Le  J 

and 


2  it  D„  2  n  Dp 

.  =  —  -  sino  .  ,  <!> .  =  ----  sino. 

d2  Xo  d  '2  Xo  1 


(153) 


Thus,  for  D9<0.05a  the  phase  shifts  <j>.  and  <j>.  and  hence  the  SINR 

C.  0  Qo  1  r\ 

are  smal 1 .  L 


The  fact  that  Ap=l  also  means  that  the  eigenvalue  spread  S  is 
large  for  all  o..  For  Dp=0.05AQ,  we  have 


S  =  s 
min  max 


2000. 


As  the  separation  increases  from  0.1aq  to  0.3aq  in  Figure  13, 
we  see  that  Ap  also  increases.  In  other  words,  we  have  better  SINR 
as  well  as  decreased  eigenvalue  spread  with  increased  element  separa¬ 
tion.  From  the  figure,  we  see  that  with  a  spacing  0.2ao,  the  two- 
element  array  is  able  to  null  interference  adequately  in  the  sector 
180°<o .<360°  for  0^45°  if  the  minimum  required  output  SINR  from  the 
array  is  not  higher  than  0.4  (-4  dB). 

The  next  graph.  Figure  14,  shows  Ap  for  the  same  conditions  as 
in  Figure  13  except  that  Dp=0.4\Q,  0.5aq  and  0.6ao.  The  SINR  improves 
in  the  sector  0">n.-180"  with  increased  separation  but  degrades  in  the 
sector  180°-'(k  -360°.  In  particular,  for  Dp=0.6\o,  the  SINR  is  low 
when  the  interference  is  coming  from  240°- o .  300°.  for  still  larger 
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Figure  14.  Effect  of  element  spacing  (0^0.4^,  0.5,*o  and  0.6aq). 

SNA-1.,  I NR- 1000. ,  BD^O.OOO,  BI-0.000, 

THETA  D-45.  DEGREES;  2-ELEMENT  ARRAY. 


separations,  the  array  starts  to  have  grating  nulls,  figure  15 
illustrates  this  phenomenon.  The  graph  shows  a^  for  Dp=0.7 \q,0.8' 
and  0.9aq,  respectively.  For  D^-0.7ao,  there  are  two  grating  nulls 
around  0^=225°  and  o.=3l5".  As  the  separation  is  increased,  the 
grating  nulls  move  in  opposite  directions  as  indicated.  In  Figure  16, 
*2  is  plotted  with  0^=1  •  anc*  1-7aq.  e*hibits  fast  fluctuations 
between  its  extrema  and  the  SINR  has  sharper  lobes  and  more  grating 
nulls. 

In  conclusion,  for  two  isotropic  elements,  the  number  of  grating 
nulls  and  the  number  of  minima  and/or  maxima  increase  with  increased 
element  separation.  With  very  small  separation,  such  as  lor  U,,  n.;v>  (), 
the  eigenvalue  spread  is  roughly  Nr.  j  which  is  high  for  strong  inter¬ 
ference  When  the  separation  is  large  enough  to  produce  grating 
nulls,  it  also  gives  large  eigenvalue  spread. 


64 


D2  =  0.9Xo 


1  p  1  1  1  1  | 

90.  180.  270.  36 

THETR  I  (DEG.) 

Figure  15.  Effect  of  element  spacing  (D?=0.7A  , 
0.8XQ  and  0.9*Q). 

SNR=1 . ,  INR= 1 000 . ,  BD=0.000, 
BI=0.000,  THETA  D=45.  DEGREES, 

2 -ELEMENT  ARRAY. 


We  now  show  how  varies  with  interelement  spacing  for  a  three 
element  array.  In  order  to  make  comparisons  among  various  cases,  we 
first  fix  the  total  lenqth  of  the  array,  D^,  and  then  change  by 
moving  the  center  element. 

Figure  17  shows  for  0^=0. 1aq,  0.3aq  and  0.5Aq  with  kept 
constant  at  a  .  From  the  figure,  we  see  the  least  desirable  choice 

is  Do=0.1>„  because  it  gives  very  low  SINR  around  o.=220°  and  o.=340°. 

2  o  J  i  i 

Considering  the  SINR  output  in  the  sector  00<_o.<1300,  the  three  cases 
give  nearly  identical  performance.  However,  in  the  sector  180°<o.< 
360°,  the  case  with  D,,=0.3ao  does  not  have  as  deep  as  a  null  as  the 
case  D2=0.5aq  or  D^O.IAq.  This  fact  suggests  that  the  equally  spaced 
array  may  not  be  the  best  choice  in  certain  conditions. 

Next,  we  keep  D2  fixed  and  vary  the  total  length  D^.  We  arbi¬ 
trarily  choose  02=0. 3Aq  and  let  03  change  from  0.5aq  to  1 . 2> 0  and  to 
2aq.  The  eigenvalue,  a2,  for  this  three  cases  is  shown  in  Figure  18. 
In  general,  increasing  the  total  spacing  gives  higher  SINR  and  gives 
sharper  nulls  around  o_.=o^  and  o.  =  180°-o^.  However,  from  Figure  18 
we  see  that  around  points  A  and  A'  the  case  D  =2a  gives  lower  SINR 
than  the  case  with  D2=1.2aq.  Thus,  with  increased  total  spacing  the 
array  will  have  additional  SINR  drops  such  as  those  around  A  and  A’ 
in  the  fiqure. 

0 •  The  Effect  of  Element  Patterns 

In  the  previous  sections,  we  have  assumed  that  the  elements  are 
isotropic.  We  shall  now  illustrate  how  element  patterns  affect  the 
eigenvalues. 

With  isotropic  elements,  and  with  the  interference  much  stronger 
than  the  desired  signal,  the  largest  eigenvalue,  \j,  exhibits  little 
variation  as  the  signal  angles  vary.  For  non-isotropic  elements, 
however,  a^  is  no  longer  constant.  Also,  the  behavior  of  the  middle 
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Figure  18.  Effect  of  element  spacings  (D2=0.3xo,  03  is  varied) 

SNR=1 . ,  INR=1000. ,  BD=0. 000,  BI=0.000 

THETA  0=45.  DEGREES 


3-ELEMENT  ARRAY. 
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eigenvalue,  x^,  is  changed.  He  shall  use  a  three  element  array  of 
dipoles  spaced  a  half-wavelength  apart  to  illustrate  these  effects. 

Let  us  assume  each  of  the  elements  is  a  short  dipole  with  a 
cosine  pattern  of  the  form 

fjo)  =  cos(o+*p) 

for  «=1,2,3.  Initially,  we  shall  assume  all  three  elements  have 
their  pattern  maxima  at  broadside,  i.e.. 


for  £=1,2,3.  Then,  all  elements  also  have  nulls  at  0=90°  and  270°. 
First  we  consider  the  case  where  the  interference  is  weak,  i.e.,  let 
Cd=l  and  =  10.  Also,  we  let  0^=45°  as  usual.  The  resultant  eigen¬ 
values  (X^jXg  and  are  shown  in  Figure  19. 

Comparing  Figure  19  with  Figure  4  (calculated  for  isotropic 
elements),  we  see  that  the  first  difference  is  in  the  largest  eigen¬ 
value,  X-| .  The  range  of  variation  of  in  Figure  19  is  much  larger 
than  that  in  Figure  4.  With  cosine  elements,  we  have 


X,  =31,  x  =  X  (o.=90°,  or  270°)  =  2.5  . 

1  T  ' 


Notice  that  x.  occurs  when  the  interference  is  arriving  in  the 
min 

pattern  nulls  at  n.=90"  or  270°.  In  the  vicinity  of  these  two  inter¬ 
ference  angles,  we  have 


T  * 

iuTui 


so  the  array  behaves  as  if  only  one  (the  desired)  signal  is  present. 

?  3  5 

Hence  x^  is  just  1  +  3  (j|f(°cj)|  =1  +  ^  =  2  checks  with  the  figure. 
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Figure  19.  Eigenvalues  vs.  o..  for  dipole  (cosine)  elements 

SRI=1 .0 
INR=10. 

8D=0.0 

BI=0.0 

THETA  D=45  DEGREES 

3-ELEMENT  ARRAY  WITH  HALF-WAVELENGTH  SPACINGS. 
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1 


From  Equation  (98),  we  know  that  for  isotropic  elements,  A-)  . 
occurs  when  the  two  signal  vectors  are  orthogonal,  i.e.,  when 


It  is  evident  that  this  condition  is  satisfied  for  o.=90°  or  <>.=270" 
with  cosine  elements  (although  we  might  not  describe  this  relation 
as  orthogonality). 

Considering  next  the  second  eigenvalue,  we  see  that  ^  has 
additional  minima  around  o.=90°  and  270°.  These  minima  occur  when¬ 
ever  the  array  behaves  as  if  only  one  signal  is  present.  We  have  from 
Figure  19  that 

A2  |  A2  .  =  1 

max  min 

Theoretically,  from  Equation  (60)  and  (62)  we  know  that  when 

T  *  T  *  5 

UjUd  :  0  ^  will  reach  its  maximum  which  is  ^ 

for  it  =45°. 
d 

Note  that  the  smallest  eigenvalue,  *3,  remains  at  unity.  This 
eigenvalue  is  not  affected  by  the  change  in  element  patterns. 

The  previous  example  shows  a^ ,  A^  and  A^  for  cosine  elements 

with  weak  interference  (o - =1 0) .  We  now  show  how  A^  and  A^  behave  for 

cosine  elements  with  strong  interference.  Suppose  we  now  set  f ..  =  1 000 

with  other  parameters  the  same  as  in  Figure  19.  We  first  show  a^  in 

Figure  20.  From  the  figure,  we  see  that  with  cosine  elements  ^  has 

a  range  of  variation  far  larger  than  that  for  isotropic  elements. 

Because  of  the  pattern  nulls  around  <<.=90"  and  270",  A  varies  between 

3  1 

a  high  value  of  approximately  3x10  to  a  low  of  2.5.  Thus,  is  very 
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Figure  20.  A-|  vs  0j  for  isotropic  elements  and  dipole 

(cosine)  elements. 
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sensitive  to  the  element  patterns  and,  of  course,  can  no  longer  be 
viewed  as  a  constant. 

Figure  21  shows  ^  for  both  isotropic  and  cosine  pattern  ele¬ 
ments.  In  general,  ^  has  the  same  shape  for  both  cases  except 
around  the  pattern  nulls  at  tK=90°  and  270".  Within  the  region 
around  these  pattern  nulls  (indicated  by  A  and  A'),  the  interference 
is  nulled  by  the  element  patterns.  I.e.,  the  array  behaves  as  if 
only  one  (the  desired)  signal  is  present. 

Notice  that  around  these  null  regions  (A  and  A1),  A„-l  is  not 

T  *  c  T  * 

related  to  the  output  SINR  because  the  condition  £..U.U.j  '^d^d  1S 

not  met. 

Next,  consider  what  happens  if  the  element  patterns  are  rotated 
so  the  elements  have  beam  maxima  and  nulls  in  different  directions. 
Rotating  the  elements  eliminates  the  symmetry  null  at  0^=180"-'^. 

For  example,  in  the  element  patterns  f|t(o)=cos(n+i^) ,  let  us  set 

<^=-60°,  -2=0o  and  a3=60°. 

a.j  for  this  situation  is  shown  in  Figure  22.  Notice  that  the  range 
of  variation  is  smaller  than  that  for  the  unrotated  pattern  case. 

We  have  from  the  figure 

A  =  1502.5  and  a  =  1501. 
max  min 

A ■)  .  does  not  drop  to  a  low  value  because  when  the  interference  is 
'min  v 

nulled  by  one  element,  it  can  still  be  picked  up  by  the  other  two 
elements  since  the  elements  have  different  null  directions.  is 
essentially  constant  in  this  case. 

In  Figure  23,  we  show  A^  for  the  rotated  elements.  Notice  that 
the  only  angle  where  ^  reaches  its  minimum  is  e .  =  i^-45'' .  For  all 
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Figure  21.  Effect  of  element  patterns. 
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Figure  22  vs  o.-  with  rotated  elements. 
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Figure  23.  A£  vs.  Oj  with  rotated  elements. 
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other  values  of  (K  we  have  A^^l.  From  Equation  (152),  we  know  that 
the  array  output  SINR  has  only  one  substantial  dip  at  o.=o^.  The 
corresponding  dip  due  to  symmetry  at  o.=180°-«j  has  been  eliminated 
by  the  rotation  of  the  elements. 


Because  a  is  independent  of  the  element  patterns,  we  know  that 
the  smallest  eigenvalue  remains  at  unity. 
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t .  The  Effect  of  Signal  Bandwidth 

The  previous  sections  have  discussed  the  eigenvalues  for  CW 
signals.  We  shall  now  study  the  effect  of  signal  bandwidth  on  the 
eigenvalues. 


In  Equation  (43),  the  element  of  the  covariance  matrix  ■ 

was  given  for  the  non-zero  bandwidth  case.  We  first  normalize  this 

2 

equation  with  respect  to  the  noise  power  ,  i.e., 

n  n  j(V"fd> 

■:'im  *  fl("d)V"d>''d  s1nc  U  Bd(V*d  ’J  e 

”  JO,  ) 


+  f  „  (o .  )f  (o sine 
i '  m  l  i 


[?  Bi<*u-*i>] 

w  m  J 


lc  1 

m 


+  >s 


k’m 


From  these  matrix  elements,  we  have  obtained  the  normalized  eigen¬ 
values  of  ‘t*’  by  a  subroutine  in  the  IBM  Scientific  Subroutine  Package 
[9],  The  subroutine  is  based  on  the  Jacobi  method  (modified  by  Von 
Neumann[10,l 1 ] )  We  shall  first  show  the  effect  of  desired  signal 
bandwidth  on  the  eigenvalues  and  then  the  effect  of  interference 
bandwidth.  We  consider  again  an  array  with  three  isotropic  elements 
a  half  wavelength  apart.  The  signal-to-noise  ratios  are  assumed  to 

be  £,.=1  and  ;',-=1000  with  o,=45n. 
d  l  d 

We  first  show  how  Ap  varies  for  different  desired  signal  band- 

widths,  Bj,  in  Figure  24.  B^  has  been  increased  from  zero  to  one  in 

steps  of  0.2  and  is  kept  at  zero.  In  the  sector  25°  <y  155°,  we 

see  that  increases  with  increased  B.,  so  the  time  constant  of  the 
2  d 

array  associated  with  'p  becomes  smaller.  On  the  other  hand,  for 
interference  coining  from  outside  that  sector,  ^  decreases  with  in¬ 
creased  B.  and  the  associated  time  constant  is  larger.  The  combined 
d 

effect  of  the  above  two  results  is  that  the  range  of  variation  of  'p 
is  reduced  with  increased  B^.  For  example,  with  B^  0.4,  we  have 
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^  =  3.88  and  a  -1.12 

max  min 

so  the  range  of  a?  has  been  reduced  to  2.76  (an  87  change  from  3)  with 

407  change  in  the  desired  signal  bandwidth.  Notice  that  a?  .  is  no 

6mi  n 

longer  unity  for  non-zero  Bj. 

We  now  consider  With  non-zero  bandwidth,  a  is  no  longer 
constant.  In  Figure  25,  we  show  how  a^  increases  with  increased  B^. 
With  Bd=0.4,  we  find 

A  =  1.15  A  =  1 .05  . 

max  min 

The  time  constant  associated  with  a^  decreases  with  increased  B^. 

Because  the  sum  of  the  three  eigenvalues  is  a  constant,  we  can 

obtain  the  behavior  of  A^  with  B^  from  Figures  24  and  25.  However, 
since  the  changes  in  both  A^  and  A^  (hence  in  are  very  small  (  5 
for  Bd=l)  compared  to  the  magnitude  of  a^ ,  which  is  around  3000,  it 

turns  out  that  A^  is  still  essentially  constant  for  B^  1 . 

Now  let  us  discuss  the  effect  of  interference  bandwidth,  . 

In  Figure  26,  we  have  calculated  a^  for  various  B.. .  It  is  seen 
that  A^  is  more  sensitive  to  Bi  than  to  B^  because  a  small  (57) 
increase  in  B^  causes  a^  to  increase  significantly  while  a  large 
(40%)  increase  in  B^  makes  a^  change  only  slightly  (87).  This 
difference  is  due  to  the  fact  that  the  interference  is  very  much 
stronger  than  the  desired  signal. 

From  the  figure,  we  see  that  ^  increases  with  B..  Also,  notice 
that  when  the  interference  is  coming  from  broadside  ('  •"O'  or  180"), 

Aj,  has  a  low  value.  This  phenomenon  occurs  because  for  these  arrival 
angles  there  is  no  interelement  delay  for  the  interference  so  the 
interference  has  the  same  effect  on  the  covariance  matrix  as  a  CW 
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Figure  25.  Effect  of  Bd  on  x^. 

SNR=1 . 

INR=1000. 

THETA  D=45. 

B I =0 . 00 
BD  is  changed. 


signal.  From  the  previous  result  in  Equation  (90)  we  know  that  a 

T  *  ^ 

is  bounded  by  for  CW  signals,  so  a is  approximately  4  at 

0^=0°  and  180°,  as  shown 


The  effect  of  B^  on  the  small  eigenvalue,  a^,  is  shown  in  Figure 
27.  As  increases,  a^  also  increases  and  exhibits  more  complicated 
behavior. 
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SECOND  EIGENVALUE  ,  X2 
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Figure  26  Effect  of  on 

SNR=1 

INR=1000 

Rd-0 

Bi  is  varied. 
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Figure  27.  Effect  of  B.  On  A.. 

SNR=1 . 

INR=1000. 

BD=0. 00 
THETA  D=45.0 
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Since  and  increase  with  and  the  sum  of  the  three  eigen¬ 
values  is  a  constant,  we  conclude  that  x^  must  decrease  with  increased 

B-.  However,  the  percentage  change  in  A,  caused  by  B  is  still  small 
1  1  i 

because  x^  is  large.  For  example,  with  B-=0.2,  we  see  from  Figures  24 

and  25  at  o.=90° 
i 

X£  =  66  ,  x_  =  1 .2, 

so  from  x-j+*2+x3  =  3006  we  know 
A,  =  2938.8. 

Thus  X  decreases  from  its  value  of  3002.6  for  CW  to  2938.8  with  20 
percent  interference  bandwidth.  The  percentage  change  is  less  than 
2  percent,  i.e., 

3002.6-2938.6  <  n 
3002.6 

Therefore,  the  effect  of  B^  on  X^  is  rather  small. 

In  conclusion,  increased  interference  bandwidth  causes  x^  and 
A^  to  increase  and  x.|  to  decrease.  The  time  constants  associated  with 
x2  and  x3  will  decrease  with  B.  and  that  associated  with  x^  will  in¬ 
crease.  Desired  signal  bandwidth  causes  the  range  of  variation  of  x^ 
to  decrease  and  that  of  x^  to  increase.  In  addition,  x^  is  little 
affected  by  both  Bd  and  Bi  provided  that  Bd<l  and  B^O.2. 
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CHAPTER  IV 
CONCLUSIONS 


In  this  report,  we  have  first  derived  the  eigenvalues  of  the 
covariance  matrix  for  an  N  element  adaptive  array  with  a  CW  desired 
signal  and  a  CW  interference.  There  are  N-2  constant  eigenvalues  in 
this  case.  The  remaining  two  eigenvalues  depend  on  the  signal  envir¬ 
onment  and  have  been  shown  to  vary  within  certain  bounds.  Hence  the 
eigenvalue  spread,  defined  as  the  ratio  between  the  two  eigenvalues, 
is  also  bounded,  as  in  Equation  (110). 

Furthermore,  we  have  discussed  the  effects  of  various  signal 
and  array  parameters  on  the  eigenvalues.  It  has  been  shown  that  when 
the  interference  is  strong,  the  largest  eigenvalue  is  essentially 
constant  for  isotropic  elements.  This  large  eigenvalue  exhibits  little 
percentage  change  as  the  signal  angles  are  varied.  The  large  eigen¬ 
value  remains  constant  regardless  of  element  spacing  or  signal  band¬ 
width.  It  does  depend  on  the  input  INR,  the  number  of  elements  and 
the  element  patterns,  however. 

We  have  also  shown  that  the  second  eigenvalue,  x^,  approaches 
a  limit  when  r,.  is  increased.  For  large  £.  and  CW  signals,  the  array 
output  SINR  is  one  less  than  this  limit,  as  shown  in  Equation  (152). 

The  higher  r..  the  better  the  approximation.  As  the  input  SNR,  r,d,  is 
increased,  the  range  of  variation  of  x^  becomes  larger,  x^  verV 
sensitive  to  element  spacings  and  element  patterns.  By  manipulating 
element  spacings  and  patterns,  we  can  modify  the  behavior  of  x^  and 
also  the  array  output  SINR.  In  addition,  x^  depends  very  much  on  B^; 
a  small  increase  in  B.  causes  x^  to  increase  significantly,  as  demon¬ 
strated  in  Figure  26.  In  the  case  of  greatest  interest  where  is 
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much  lower  than  the  effect  of  desired  signal  bandwidth  on  is 
much  less. 

The  remaining  N-2  eigenvalues  will  differ  from  unity  only  when 
non-zero  bandwidth  signals  are  present.  Also,  the  effect  of  on  x^ 
is  larger  than  that  of  when  the  interference  is  much  stronger  than 
the  desired  signal . 
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